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Abstract 


The  circadian  system  or  "biological  clock"  orchestrates  the  temporal  organization  of 
several  physiological  and  metabolic  processes  over  a  24-hour  period  in  most  living 
organisms.  In  mammals,  the  circadian  timing  system  has  a  hierarchical  organization 
composed  of  the  central  pacemaker  in  the  suprachiasmatic  nucleus  (SCN)  which 
coordinates  the  oscillating  activity  of  peripheral  clocks  that  are  present  in  almost  all 
tissues.  At  the  molecular  level,  circadian  oscillations  in  both  SCN  and  peripheral  cells  are 
based  on  negative  autoregulation  of  gene  expression  involving  essentially  the  same  core 
clock  components.  At  the  tissue  level,  the  SCN  consists  of  a  coupled  network  of 
heterogeneous  neurons  that  fire  synchronously  through  neurotransmitter  coupling  and  are 
entrained  to  24-hour  period  through  daily  light-dark  cycles.  At  the  organismal  level,  the 
SCN  synchronizes  peripheral  tissue  clocks  via  hormonal  and  indirect  cues  such  as 
feeding  time,  probably  through  modifications  of  metabolite  availability.  More  recently, 
experiments  have  shown  that  interactions  between  circadian  and  metabolic  systems  are 
bidirectional  in  that  metabolites  can  feedback  into  the  clock.  These  findings  emphasize 
the  need  of  a  multiscale  integrative  approach  to  understanding  circadian  regulatory 
processes.  Here,  we  adopt  a  quantitative  systems  biology  approach  that  uncovers  critical 
properties  of  the  circadian  oscillator  at  various  scales  (neurons,  SCN,  periphery). 

Keywords:  computational  and  theoretical  biology,  biological  rhythms,  multiscale 
modeling,  metabolism,  synchrony  and  entrainment 


Scientific  Progress  and  Accomplishments 

The  circadian  system  or  "biological  clock"  orchestrates  the  temporal  organization  of 
several  physiological  and  metabolic  processes  over  a  24-hour  period  in  most  living 
organisms  from  bacteria  to  plants  and  mammals.  This  highly  conserved  system  enables 
the  organism  to  adapt  to  common  daily  and  seasonal  changes,  such  as  the  day-night  cycle 
and  food  availability.  In  mammals,  the  circadian  timing  system  is  composed  of  the 
central  pacemaker,  which  resides  within  the  suprachiasmatic  nucleus  (SCN)  and 
coordinates  the  oscillating  activity  of  peripheral  clocks  that  are  present  in  almost  all 
tissues.  At  the  molecular  level,  circadian  oscillations  in  both  SCN  and  peripheral  cells  are 
based  on  negative  autoregulation  of  gene  expression  involving  essentially  the  same  core 
clock  components.  At  the  tissue  level,  the  SCN  consists  of  a  coupled  network  of 
heterogeneous  neurons  that  fire  synchronously  through  neurotransmitter  coupling  and  are 
entrained  to  24-hour  period  through  daily  light-dark  cycles.  At  the  organismal  level,  the 
SCN  synchronizes  peripheral  tissue  clocks  via  humoral  signals  and  indirect  cues  such  as 
feeding  time,  probably  through  modifications  of  nutrient  availability.  More  recently, 
experiments  have  shown  that  interactions  between  circadian  and  metabolic  systems  are 
bidirectional  in  that  the  clock  not  only  influences  metabolic  processes  but  also  metabolic 
signals  feedback  into  the  circadian  machinery.  Taken  together,  the  circadian  system  has  a 
hierarchical,  multiscale  organization  and  integrative  quantitative  approaches  are  crucial  to 
understanding  circadian  regulatory  processes.  Here,  the  opportunities  of  a  multiscale 
integrative  approach  will  be  discussed  for  systems-level  analysis  and  modeling  in 
circadian  physiology.  The  unique  aspects  of  this  work,  which  are  discussed  in  great  detail 
in  the  relevant  publications  [1-9],  are  summarized  below  and  include: 

(i)  development  of  a  cell  autonomous  (single  cell)  mathematical  model  of  a 
circadian  neuron  which  resolves  the  experimentally  observed  discrepancies 
between  the  tissue  scale  and  the  cellular  scale; 

(ii)  generation  of  a  tissue  level  model  by  coupling  circadian  neurons  using  the 
vasoactive  intestinal  polypeptide  (VIP)  predicting  tissue-level  phenotypes; 

(iii)  development  of  nonlinear  model  predictive  control  strategies  to  control 
circadian  behavior; 

(iv)  development  of  novel  sensitivity  and  model  reduction  techniques  used  to 
assess  both  the  robustness  of  the  model  to  parameter  perturbation  and 
reduce  the  dimensionality  scale  for  computational  feasibility; 

(v)  development  of  “integrate-and-fire”  models  to  represent  synchronized 
population-scale  phenomena  in  many  other  biological  systems  including  the 
annual  spawning  of  certain  species  of  corals  and; 

(vi)  derivation  of  rigorous  analytical  metrics  to  understand  general  principles  of 
biological  synchrony  and  circadian  entrainment. 

Our  most  recent  activities  continue  to  probe  at  the  sources  of  regulation  that  give  rise  to 
robust  oscillations  and  highly  precise  periods  in  the  mammalian  “biological”  clock. 


Firstly,  upon  recognition  of  the  reciprocal  interaction  between  the  circadian  clock  and 
metabolism,  we  have  developed  a  novel  circadian-enzymatic  model  that  quantifies  the 
metabolic  regulation  of  the  core  circadian  oscillator.  Such  model  generates  new 
experimentally  testable  hypotheses  that  can  explain  a  contradiction  in  the  metabolic/clock 
connection.  Secondly,  we  have  also  begun  to  obtain  theoretical  results  in  the  regulation  of 
the  collective  period  in  a  coupled  circadian  oscillatory  network.  Using  an  autorepression- 
based  gene  regulatory  network,  such  analytical  insights  indicate  that  in  the  presence  of  a 
coupling  delay,  the  strength  of  the  coupling  (intercellular)  neurotransmitter  regulates  the 
collective  period.  Further  details  of  our  previous  and  ongoing  research  activities  are 
discussed  below. 

Modeling  cell-autonomous  oscillations  of  a  circadian  neuron 

The  fundamental  structure  of  the  biochemical  circadian  oscillator  in  mice  has  been 
known  for  roughly  2-3  decades.  This  has  led  to  the  construction  of  various  mathematical 
models  of  the  mammalian  clock,  both  deterministic  and  stochastic.  These  models  were 
constructed  from  tissue-level  data  and  therefore  describe  the  behavior  of  the  molecular 
clock  in  communicating  cells.  Recent  experiments,  however,  show  that  the  gross 
phenotype  of  the  system  (i.e.,  rhythmicity  vs.  arrhythmicity)  is  often  different  in  gene 
knockouts  between  autonomous  and  non-autonomous  cells  [10]. 

To  maintain  proper  functionality, 
individual  clock  components  must  rise  and 
fall  in  fixed  phase  relationships  with  one 
another,  alternately  turning  on  and  off  the 
genes  in  the  network.  We  built  a  network 
structure  [5],  Figure  1,  from  the  most 
recent  experimental  data,  incorporating 
eight  key  genes  -  more  than  are  found  in 
any  previous  mathematical  model  of  the 
mouse  clock.  We  then  used  the  known 
phase  relationships  among  mRNAs  and 
between  mRNAs  and  proteins,  Figure  2,  to 
construct  a  cost  function  and  found  a  best- 
fit  parameter  set  for  the  model  by 
minimizing  this  cost  through  an  iterative 
evolutionary  strategy.  The  resulting  model  was  validated  against  the  gross  phenotypes  of 
seven  gene  knockouts  and  two  double  gene  knockouts.  For  all  nine  cases,  the  model 
correctly  predicts  the  phenotype.  Additional  validation  was  achieved  via  comparison  to 
the  known  molecular  biology  of  the  system. 
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Figure  1:  Network  topology  of  a  cell-autonomous 
circadian  neuron.  Interacting  components  involve  core 
clock  feedback  loops  (e.g.  PER/CRY)  and  secondary 
(auxiliary)  loops  (e.g.  ROR/REV-RB)  regulating  the 
circadian  clock. 


Further,  we  assessed  the  robustness  of  our 
wild-type  model  to  maintain  periodicity  in 
response  to  parameter  perturbation,  a 
requirement  to  ensure  that  our 
parameterization  is  not  so  restrictive  as  to 
be  unrealistic.  We  found  that  the  model 
behavior  is  maintained  in  98%  of 
alternative  models.  Taken  these  together, 
although  this  model  was  calibrated  using 
only  the  phase  relationships  among  mRNA 
and  protein  concentrations,  it  accurately 
predicts  all  the  single  cell  phenotypes  of 
gene  knockout  mutations  for  which 
experimental  data  exist. 

Tissue-level  modeling  of  coupled  autonomous  circadian  neurons 

While  the  gene -regulatory  feedback  networks  in  each  SCN  pacemaker  produce  "noisy" 
self-sustained  rhythms,  robust  approximately  24-hour  output  is  generated  from  the  SCN 
through  the  coupling  of  thousands  of  pacemakers  through  neurotransmitter  coupling  and 
is  entrained  to  24  hour  period  through  daily  light-dark  cycles. 

Synchronization  among  neurons  in  the  SCN 
is  achieved,  at  least  in  part,  through  circadian 
release  of  the  VIP  and  its  reception  by  the 
VPAC2  receptor,  leading  to  the  initiation  of  a 
signal  cascade  involving  activation  of  G- 
protein,  internalization  of  the  VPAC2 
receptor,  activation  of  adenylate  cyclase, 
production  of  cyclic  AMP  (cAMP), 
activation  of  protein  kinase  A  (PKA), 
phosphorylation  of  cAMP-response  element 
binding  protein  (CREB),  and  finally,  up- 
regulation  of  one  or  more  Per  genes. 

Circadian  oscillations  have  been  observed  in 
cAMP  and  calcium  levels  and  also  in  CREB- 
mediated  gene  expression  in  the  SCN.  Light, 
the  principal  agent  responsible  for 
entrainment  of  the  internal  clock  to  the  environment,  is  transduced  on  the  retina  and 
channeled  to  the  SCN  via  the  retinohypothalamic  tract,  where  its  effects  ultimately  lead 
to  up-regulation  of  one  or  more  period  (Per)  genes.  The  path  that  leads  to  light-induced 
up-regulation  may  be  similar  or  identical  to  that  which  synchronizes  the  SCN  neurons  to 
one  another,  since  VIP  release  appears  to  be  enhanced  by  light. 


Figure  3:  Schematic  of  the  coupled  (tissue-level) 
model  of  the  mammalian  circadian  clock. 

Communication  among  cells  is  achieved  via  circadian 
release  of  the  vasoactive  intestinal  polypeptide  (VIP), 
which  is  received  and  transduced  by  the  VPAC2 
receptor.  Figure  is  adapted  from  [11]. 
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Figure  2:  Dynamic  profiles  of  the  mammalian 
circadian  clock.  Blue  and  red  dots  represent 
experimentally  observed  phase  relationship  between 
mRNA  and  proteins  respectively. 


We  incorporated  the  signaling  cascade 
described  above  and  shown  in  Figure  3  using 
the  mechanistic  representation  used  by  To  et 
al.  [12]  to  couple  a  5x5  grid  of  our  cell 
autonomous  model  in  Mirsky  et  al.  [5].  We 
found  parameters  for  the  22  parameters  in  the 
VIP  cascade  that  produced  circadian 
oscillations  in  the  tissue-level  model  while 
retaining  all  the  original  parameters  for  the  cell 
autonomous  model.  These  oscillations  (see 
Figure  4)  are  produced  by  synchronized 
oscillations  in  the  25  cells  with  a  period  of  30 
hours  and  moreover,  these  oscillations  are 
robust  in  the  presence  of  parameter 
heterogeneity  in  the  coupled  cells.  As 
mentioned  earlier,  the  cell  autonomous  model 
can  predict  all  the  single-cell  knockout 
phenotypes  observed  in  experiments. 
However,  some  knockout  phenotypes  at  the  tissue-level  are  different  from  the  single-cell 
phenotypes  [10].  Unfortunately,  the  deterministic  model  showed  lack  of  oscillations  at 
the  current  parameter  values  with  Perl,  Per2  and  Cry  I ,  which  give  robust  rhythms  in  the 
tissue.  Nevertheless,  we  tested  the  knockdown  (instead  of  knockout)  of  these  3  genes  and 
were  still  able  to  obtain  rhythmic  behavior  as  seen  experimentally. 

Synchrony  and  entrainment  properties  of  robust  biochemical  oscillators 

In  addition  to  mathematical  modeling  (see 
above),  we  also  employed  advanced  systems 
theoretic  tools  (e.g.  control  systems  analysis 
algorithms)  in  order  to  better  understand, 
characterize  and  control  the  circadian 
timekeeping  system  [2,  3].  Precise 

synchronization  of  these  rhythms  is  an 
essential  part  of  circadian  organization, 
which  is  hierarchically  organized  into  the 
following  steps:  (i)  periodic  light  detected  by 
the  eyes  entrains  (synchronizes)  -20,000 
neurons  in  the  SCN  via  neuronal  signals;  (ii) 
within  SCN  tissue,  individual  neurons 
synchronize  each  other  via  VIP  coupling 
and/or  gap  junctions  and  subsequently  (iii) 
the  SCN  synchronizes  the  other  peripheral 
tissue  clocks  via  humoral/indirect  signals  and 
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Figure  5:  State  dynamics  and  light  profiles  used 
for  circadian  entrainment.  (Top)  State  dynamics  as 
they  converge  to  the  reference  trajectory  (bold  solid 
line);  (Bottom)  Light  profiles  of  duration  1  (dotted 
lines),  2  (dot-dashed  lines)  and  3  hours  (dashed  lines) 
used  to  reset  the  phase  differences;  these  sequences 
also  converge  to  nominal  light/dark  cycles  (bold  solid 
line). 
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Figure  4:  Time  courses  for  a  5  x  5  grid  of 
coupled  cell-autonomous  circadian  clocks.  Top 

plots  show  Perl  mRNA  concentrations  for  all  25 
cells,  middle  plots  show  average  Perl  mRNA 
concentration  across  all  25  cells,  and  bottom  plots 
show  synchronization  indices  (Sis). 


such  interaction  is  reciprocal  in  nature.  Since  little  is  known  about  this  communication, 
our  study  investigated  the  coupling  mechanisms  that  give  rise  to  synchronized  networks 
in  the  SCN.  Using  deterministic  and  discrete  stochastic  models  in  both  the  single-cell  and 
network  setting,  we  investigated  optimal  control  strategies  for  light-induced  circadian 
phase  entrainment  as  illustrated  in  Figure  5.  Specifically,  we  tested  the  efficacy  of  the 
optimization  (model  predictive  control,  MPC)  algorithm  with  respect  to  various  control 
(light)  inputs  as  applied  to  a  nonlinear  mammalian  model.  Further,  our  investigation  of 
the  Drosophila  clock  model  revealed  that  only  82  (out  of  960)  cellular  coupling 
mechanisms  would  provide  synchrony  and  therefore  significantly  decreasing  the  number 
of  candidates  to  be  tested  experimentally.  Thus,  the  analysis  of  phase  characteristics  may 
facilitate  the  identification  of  relevant  coupling  mechanisms  for  the  ensemble  of  noisy 
(stochastic)  circadian  clocks. 

Recently,  our  research  efforts  in  systems 
biology  have  also  expanded  into  the 
application  domains  of  ecology,  with  the 
broad  technical  theme  of  understanding 
synchronized  population-scale  phenomena, 
such  as  coral  spawning,  using  coupled  and 
driven  oscillator  models  [1,4].  These  studies 
identify  the  sources  of  coherence  in  the 
annual  spawning  of  certain  species  of  coral, 
and  this  builds  on  the  group’s  expertise  in 
robustness  analysis  of  coupled  oscillators  in 
the  brain  that  are  responsible  for  circadian 
rhythms.  More  specifically,  we  developed  a 
bioenergetic  model  integrate-and-fire  model, 
Figure  6  that  reveals  how  annual  insolation 
rhythms  can  entrain  the  gametogenetic  cycles 
in  tropical  hermatypic  corals  to  the 
appropriate  spawning  season  since 
photosynthate  is  their  primary  source  of 
energy.  In  the  presence  of  short-term 
fluctuations  in  the  energy  input,  a  feedback 
regulatory  mechanism  is  likely  required  to 
achieve  coherence  of  spawning  times  to 
within  one  lunar  cycle,  in  order  for 
subsequent  signals  such  as  lunar  and  diurnal 
light  cycles  to  determine  the  “correct”  night 
of  spawning.  The  feedback  mechanism  can  also  provide  robustness  against  population 
heterogeneity  that  may  arise  due  to  genetic  and  environmental  effects. 
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Figure  6:  A  bioenergetic  integrate-and-fire  model 
of  coral  spawning.  (Top)  Allocation  to  the 
reproductive  reserves  Xt  of  energy  obtained  by  a  coral 
individual  from  its  photoautotropic  symbiont  derived 
from  insolation  St  with  fluctuations  nt  with  the 
remaining  energy  going  toward  growth,  somatic 
maintainance  and  other  reserves.  ^(Xt,  St)  represents 
the  fractional  allocation  toward  reproduction; 
(Bottom)  The  percentage  of  the  individuals  in  a 
colony  that  spawn  in  a  day  for  a  5-year  period  for  the 
simple  reserve  accumulation  under  A(Xt,  St)  =  St 


Novel  sensitivity  techniques  for  highly  variable  biochemical  systems 

Classical  sensitivity  and  model  reduction 
techniques  are  long-standing  systems  biology 
approaches  used  both  to  gain  insight  into  model 
subprocesses  and  reduce  the  computational  cost  of 
simulation  and  analysis.  For  systems  with  periodic 
and  stochastic  behavior  the  situation  is  somewhat 
challenging.  Recently,  we  developed  a 
distribution-based  sensitivity  analysis  method  [8] 
that  addresses  the  robustness  of  a  model  (system) 
to  extrinsic  noise  in  parameters  in  both 
deterministic  and  stochastic  regimes.  We 
particularly  applied  this  method  to  a  highly 
nonlinear  system,  the  Drosophila  biological  clock. 

As  seen  in  Figure  7  the  fragility  of  this  cell 
autonomous  model  is  concentrated  in  processes 
whose  rates  are  expressed  in  parameters  that  are 
not  clock  specific  (e.g.,  general  transcription, 
translation),  while  the  clock  is  robust  to 

perturbations  in  rates  of  reactions  that  are  found 
only  in  the  clock  network  (e.g.,  association  and 
dissociation  of  clock  proteins  that  form 
heterodimers).  Remarkably,  in  spite  of  being  a 
highly  nonlinear  system,  the  circadian  clock 
produced  proportional  larger  deviations  for 
increasing  parameter  perturbations.  Finally  we 
showed  that  the  circadian  clock  is  quite 
insensitive  to  extrinsic  noise  in  the  form  of 
fluctuations  in  the  number  of  cellular 
components  that  alter  reaction  rates  at  multiple 
genes  and  gene  products. 

In  addition  to  sensitivity  analysis,  we  also 
developed  a  novel  model  reduction  technique  for 
circadian  models  that  preserves  the  phase  response  behavior,  a  key  performance  criterion 
for  biological  clocks  [6,  9].  In  our  previous  studies  [2,  3]  we  have  already  shown  that 
phase  entrainment  can  serve  as  a  critical  control  objective  to  investigate  dynamics  both  at 
the  single  cell  level  and  at  the  population  level  using  coupled  oscillators.  The  proposed 
methodology  effectively  reduced  Forger  and  Peskin’s  61-state  mammalian  model  [13]  to 
13  states  -  a  fivefold  reduction  in  size  as  illustrated  in  Figure  8.  Specifically,  we 
discovered  that  four  feedback  loops  are  unnecessary  (redundant)  with  respect  to  the  phase 
response  behavior  of  the  single  cell  oscillator.  Using  a  coupled  multioscillator  circadian 
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Figure  8:  Oscillator  model  reduction 
preserving  the  phase  response  behavior.  Cost 
and  model  size  of  the  10  fittest  individuals  from 
each  22  generations  created  using  a  genetic 
algorithm  hybridized  with  a  local  deterministic 
search.  The  earliest  generations  are  shown  in 
blues  and  the  latest  are  shown  in  reds. 
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Figure  7:  Box  plots  of  sensitivities  classed 
by  parameter  type.  Global  parameters  are 
those  that  are  not  clock-specific  (e.g.,  general 
transcription,  translation),  local  parameters  are 
clock  specific  (e.g.,  association  and 
dissociation  of  clock  components),  and  mixed 
parameters  have  characteristics  of  the  other 
two  classes.  Note  that  the  clock  is  principally 
sensitive  to  perturbations  in  global  parameters. 


model,  we  further  demonstrated  that  by  preserving  the  phase  response  properties  of  a 
single  oscillator,  the  timing  behavior  at  the  population  level  is  also  preserved. 

Modeling  the  reciprocal  interaction  between  circadian  timekeeping  and  metabolism 

In  the  last  decade,  many  aspects  of  the  molecular  mechanism  of  circadian  rhythms  have 
been  elucidated,  thanks  to  the  advanced  progress  in  molecular  biology  experiments. 
While  it  has  been  a  dogma  for  years  that  the  circadian  clock  is  regulated  by  transcription 
feedback  loops,  emerging  studies  now  highlight  the  importance  of  another  loop,  an 
enzymatic  one.  Of  particular  interest  are  recent  groundbreaking  studies  [14,  15]  showing 
that  the  NAD+-dependent  enzyme  SIRT1  functions  as  a  histone  deacetylase  that  regulates 
core  clock  components  including  the  proteins  BMAL1  and  PER2.  While  these  studies 
show  that  the  histone  deacetylase  SIRT1  is  involved  in  regulating  the  circadian 
amplitude,  they  differ  in  the  response  of  gene  expression  to  loss  of  SIRT1  as  recently 
reviewed  in  [16].  Predicated  upon  this  apparent  contradiction,  which  we  refer  to  as 
“ SIRT1  paradox”,  this  study  aims  to  explore  such  paradoxical  effects  through  the 
development  of  a  novel  circadian-enzymatic  model,  schematized  in  Figure  9(left),  with 
the  perspective  of  generating  new  experimentally  testable  hypotheses. 


Figure  9:  (Left)  Network  topology  of  the  circadian-enzymatic  model.  Upper  case  letters  represent  the  interlocked 
transcriptional  and  enzymatic  feedback  loops  exerted  by  the  canonical  PER-CRY  loop  (A)  and  the  enzymatic  NAD 
loop  which  regulates  the  deacetylation  of  BMAL1  activator  (B)  and  PER-CRY  repressor  (C);  (Right)  Model 
recapitulates  the  experimentally  observed  differential  amplitude  response  due  to  lack  of  SIRT1.  Simulations  were 
generated  using  two  independent  parameters  sets.  Solid  lines  represent  the  wild  type  condition  while  dashed  lines 
represent  SIRT1  null  mutation.  Parameter  set  HI  was  selected  because  the  circadian  amplitude  increases  in  response  to 
SIRT1  deficiency,  while  parameter  set  H2  corresponds  to  a  situation  in  which  the  amplitude  decreases  in  the  SIRT1"  " 
mutant. 


The  performance  of  the  resulting  model  is  evaluated  through  its  potential  to  capture 
experimentally  observed  circadian  phenotypes  of  gene  knockouts  including  circadian 
amplitude  variation  due  to  loss  of  metabolic  activity  (e.g.  SIRT1'  '),  Figure  9(right). 
Interestingly,  the  model  recapitulates  the  differential  amplitude  response  due  to  lack  of 
SIRT1  as  experimentally  shown  in  [14,  15].  Specifically,  the  high-amplitude  circadian 
oscillations  as  shown  by  Nakahata  et  al.  [15]  is  simulated  for  stoichiometric  ratios 
between  the  activator  and  the  repressor  that  are  comparatively  smaller  than  those  giving 


rise  to  the  reduced  circadian  amplitude  as  reported  by  Asher  et  al.  [14].  In  this  case  ( H2 
phenotype)  such  high  ratios  are  attributed  to  higher  levels  of  the  acetylated  activator  and 
lower  levels  of  the  acetylated  repressor  complex.  This  demonstrates  the  crucial  role  a 
genetic  variation  in  the  background  levels  of  the  acetyltransferases  -  enzymes  regulating 
the  acetylation  levels  of  BMAL1  and  PER  proteins  -  might  play  in  modulating  their 
stoichiometry  and  thereby  the  amplitude  phenotype  under  null  enzymatic  mutations  (e.g. 
SIRTl'').  Interestingly,  upon  saturating  the  effect  of  acetylation  on  PER  while  reducing 
the  acetylation  of  BMAL1  the  model  predicts  increased  rather  than  reduced  amplitude  in 
the  absence  of  SIRTl.  Under  these  conditions  the  effect  of  SIRTl  deficiency  on  the 
negative  (PER)  component  becomes  negligible  and  its  effect  on  the  positive  (BMAL1) 
limb  dominates. 

Experimentally,  a  decrease  in  BMAL1  acetylation  would  be  analogous  to  a  decrease  in 
CLOCK  gene  expression  which  elicits  BMAL1  acetylation.  In  regard  to  the  saturation 
effect  of  PER  acetylation,  a  possible  scenario  would  be  to  overexpress  the  enzyme  that 
acetylates  PER  protein.  Given  that  such  an  enzyme  is  not  known  the  model  suggests  an 
alternative  scenario  that  involves  a  reduction  in  the  degradation  rate  of  PER  protein. 
Since  the  kinase  (CKI)  promotes  proteasomal  PER  degradation,  reducing  CKI  levels 
would  be  equivalent  to  reducing  the  degradation  of  PER  protein.  As  part  of  this  study  and 
in  collaboration  with  colleagues  at  UPenn,  we  aim  to  conduct  the  experimental 
verification  experiments  using  the  RNA  interference  (RNAi)  technology. 

Intercellular  coupling  strength  regulates  the  collective  period  of  delay-coupled 

biological  oscillators 

Most  biological  rhythms  are  generated  by  multiple  inter-coupled  cellular  oscillators  that 
operate  in  synchrony.  Experimental  studies  showed  that  in  the  presence  of  a  coupling 
delay,  the  collective  period  can  differ  significantly  from  the  individual  oscillation  period. 
The  phenomenon  is  intriguing  in  that  the  period  of  a  biological  clock  is  of  vital 
importance  and  is  fundamental  to  many 
biological  functions.  Autorepression-based 
gene  regulatory  oscillator  is  thought  to  be 
the  core  of  many  biological  clocks  such  as 
cell  cycles  and  the  vertebrate  segmentation 
clock,  and  it  is  also  thought  be  to  at  the 
heart  of  multiple  circadian  systems.  In  such 
gene  regulatory  oscillators,  the  protein 
directly  binds  to  the  regulatory  DNA  of  its 
own  gene  to  inhibit  transcription,  with  an 
explicit  delay  involved  in  the  synthesis  of 
mRNA  and  protein  (c.f.  Figure  10).  Based  on  delay-coupled  autorepression  oscillators, 
we  studied  analytically  the  influence  of  intercellular  coupling  strength  on  the  collective 
period.  Consistent  with  previously  published  models,  the  intercellular  coupling  is 
assumed  to  be  of  a  diffusive  type. 


Figure  10:  The  schematic  of  autoregression  based 
gene-regulatory  oscillators.  Tm  and  Tp  denote  the  time 
delay  in  transcription  and  translation,  respectively. 


We  prove  analytically  that  in  the  presence  of  a  coupling  delay,  the  strength  of 
intercellular  coupling  regulates  the  collective  period.  More  specifically,  we  prove  that 
depending  on  delay  magnitude,  the  collective  period  can  increase  or  decrease  with  the 
coupling  strength.  We  also  give  the  explicit  specific  region  in  parameter  space  in  which 
the  collective  period  increases/decreases  with  the  coupling  strength.  To  our  knowledge, 
this  is  the  first  time  that  the  collective  period  of  coupled  mechanism-based  oscillators  is 
analytically  obtained  and  is  rigorously  proven  to  be  regulated  by  the  intercellular 
coupling  strength.  Moreover,  the  collective  period  is  analytically  obtained  in  a  closed- 
form,  which  has  not  been  reported  in  the  literature.  The  results  are  given  in  terms  of 
biochemical  parameters,  and  may  give  guidance  in  synthetic  biological  oscillator  design. 
The  prevalence  of  autorepression  in  biological  oscillations  supports  the  potentially  broad 
application  to  other  biological  clocks. 

We  utilize  the  experimental  data  [17]  of  the  segmentation  clock  in  zebrafish  embryos  to 
verify  our  theoretical  predictions.  The  segmentation  clock  controls  the  somitogenesis  in 
vertebrate  development.  It  is  composed  of  many  cellular  gene  regulatory  oscillators  [18]. 
These  cellular  oscillators  are  located  within  the  posterior  mesoderm  (PSM)  and  Delta- 
Notch  signaling  synchronizes  the  cellular  network  (c.f.  Figure  11).  Treating  wild-type 
zebrafish  embryos  with  DAPT,  ay-secretase  inhibitor  that  suppresses  Notch  receptor 
function,  will  attenuate  coupling  strength  but  will  not  influence  the  time  delay  in  Delta- 
Notch  coupling  [17].  Using  our  theoretical  results,  we  can  analytically  predict  the 
collective  period  under  the  treatment  of  different  DAPT  concentrations,  as  shown  by  the 
red  dashed  line  in  Figure  ll(left).  The  experimental  measurement  of  collective  period 
under  different  DAPT  concentrations  is  borrowed  from  [17]  and  is  represented  by  green 
squares  and  error  bars  in  the  Figure  ll(right).  The  analytical  prediction  is  in  very  good 
quantitative  agreement  with  the  experimental  data. 


Figure  11:  (Left)  Schematic  of  the  presomitic  mesoderm  (PSM).  The  purple  area  indicates  waves  caused  by  the 
periodic  expression  of  the  her  gene;  (Right)  Experimental  confirmation  of  the  theoretical  results.  The  coupling 
strengths  ( k2 )  under  different  concentrations  of  DAPT  treatment  are  calculated  based  on  the  linear  relation  between 
coupling  strength  and  Notch  signaling  level  and  the  fact  that  DAPT  suppresses  Notch  signaling  following  Michaelis- 
Menten  kinetics  [17].  The  experimental  data  of  wide-type  zebrafish  embryos  treated  with  different  DAPT 
concentrations  are  represented  by  green  squares  and  error  bars.  The  analytic  prediction  is  in  very  good  quantitative 
agreement  with  the  experimental  data. 
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